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ABSTRACT 

We report a joint analysis of the Rossiter-McLaughlin (RM) effect with Subaru and the Kepler 
photometry for Kepler Object of Interest (KOI) 94 system. The system comprises four transiting 
planet candidates with orbital periods of 22.3 (KOI-94.01), 10.4 (KOI-94.02), 54.3 (KOI-94.03), and 3.7 
(KOI-94.04) days from the Kepler photometry. We performed the radial velocity (RV) measurement 
of the system with the Subaru 8.2 m telescope on August 10, 2012 (UT), covering a complete transit of 
KOI-94.01 for ~ 6.7 hours. The resulting RV variation due to the RM effect spectroscopically confirms 
that KOI-94.01 is indeed the transiting planet, and implies that its orbital axis is well aligned with 
the stellar spin axis; the projected spin-orbit angle A is estimated as — deg. This is the first 
measurement of the RM effect for a multiple transiting system. Remarkably, the archived Kepler 
lightcurve around BJD=2455211.5 (date in UT January 14/15, 2010) indicates a "double transit" 
event of KOI-94.01 and KOI-94.03, in which the two planets transit the stellar disk simultaneously. 
Moreover, the two planets partially overlap each other, and exhibit a "planet-planet eclipse" around 
the transit center. This provides a rare opportunity to put tight constraints on the configuration of 
the two transiting planets by joint analysis with our Subaru RM measurement. Indeed, we find that 
the projected mutual inclination of KOI-94.01 and KOI-94.03 is estimated to be 5 = —1.15° ± 0.55°. 
Implications for the migration model of multiple planet systems are also discussed. 
Subject headings: planets and satellites: individual (KOI-94) - planets and satellites: formation - 
stars: rotation - techniques: radial velocities - techniques: spectroscopic 



1. INTRODUCTION 

The proximity of giant planets to their host stars im- 
plies that such planets have gone through so-called plan- 
etary migrations. The angle between the stellar spin axis 
and planetary orbital axis is an impor tant probe to reveal 
the origin of plan etary migrations (jQueloz et al.l 120001: 
W inn et all 120051 ) . While migrations at earlier stages 
due to interactions between planets and proto-planctary 
disk are supp osed to result in sm aller values for the spin- 
orbit angles (|Lubow fc Ida! [20Toh , migrations caused by 
planet-planet scatterings or the Kozai cycle and subse- 
quent tidal interactions wit h the host star could pro- 
duce large spin-orbit angles ([Marzari fe W cidcnsc hillinJ 
2002tlFabrvckv fc Tremainell2007t IChatteriee et alJl2008t 



Naeasawa fc Idall201lD 



Observations of the Rossiter-McLaughlin (RM) ef- 
fect have been the major channel to measure the sky- 
projected sp in-orb it angles A (e.g.. lOhta et al.l 120051 : 
iNarita et al1l2007t lAlbrecht et al.ll2012l) . With previous 
dedicated campaigns, last few years have witnessed an 
unexpectedly large fraction of spin-orbit misaligned sys- 
tems, and even retrograde orbits with respe ct to the pro- 
jected stellar equator have been reported ([Narita et al.l 
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120091: IWinn et al.l 120091: iTriaud et all 120101) . indicating 
that dynamical processes such as planet-planet scatter- 
ings are indeed responsible for some systems with close-in 
giant planets. 

The spin-orbit angle for multiple transiting systems 
is particularly interesting. Such a multiplicity strongly 
suggests that the system has never experienced chaotic 
processes and rather results from a quiescent evolution 
history. A spin-orbit misalignment in a multiple plan- 
etary system implies that stellar obliquities can vary 
by processes whic h are unconnected to planetary mi- 



grations, (e.g., iLaTet all 1201 It iRogers et ail (20121 ) 
ISanchis-Ojeda et all (j2012l ) measured the spin-orbit re- 
lation for the Kepler-30 system by a precise modeling of 
spot-crossing events during planetary transits and found 
a good spin-orbit alignment in the system, supporting 
the notion that the system has evolved quiescently in 
the disk. 

In this letter, we focus on the KOI-94 system. This sys- 
tem, fistedjJiHie^rliest Kejjler_Oj2ject oflirterest (KOI) 
list (|Borucki et al.l 1201 It iBatalha et all 120121 ). has four 
transiting planet candidates around a relatively bright 
late F star (with the Kepler magnitude of K p = 12.205). 
The orbital period and radius R p normalized by the stel- 
lar radius R s for each candidate are summarized in Table 
[TJ Indeed KOI-94 is a remarkable system that experi- 
enced a "double transit" event due to the simultaneous 
transit of two planets on the stellar disk, and a "planet- 
planet eclipse" due to the partial overlap of the two plan- 
ets during the transiting phase. As we describe below, 
a joint analysis of our Subaru radial velocity data and 
the Kepler lightcurve provides a unique opportunity to 
tightly constrain the configuration of the multi-transiting 
planetary system. 
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TABLE 1 

Properties of Planetary Candidates in KOI-94 



Candidate 


Orbital Period (days) 


Rp / R s 


KOI-94.01 
94.02 
94.03 
94.04 


22.343000 ±0.000011 
10.423707 ± 0.000026 
54.31993 ± 0.00012 
3.743245 ±0.000031 


0.06856 ± 0.00012 
0.02544 ± 0.00012 
0.04058 ± 0.00013 
0.01045 ± 0.00019 



2. OBSERVATIONS AND DATA REDUCTION 

Due to the long orbital distance of KOI-94.01 as a tran- 
siting planet, the transit duration is considerably long 
(about 6.7 hours), which makes it very challenging to ob- 
serve a complete transit within a night by a ground-based 
telescope. A fortunate opportunity was given on 2012 
August 10 (UT), when a complete transit of KOI-94.01 
was observed by High Dispersion Spectrograph (HDS) in- 
stalled on the Subaru 8.2 m telescope. We also obtained 
out-of-transit spectra on 2012 July 1, 2, 3, 4, August 8, 
9 (UT), to determine the RV baseline and derive the RV 
semi-amplitude of the host star. In order to maximize 
the RV precision , we adopted the I2 a setup, and utilized 
the image sheer dTaiitsu et al.|[20~12T ) with the Iodine cell 
for the precise RV calibration, achieving the spectral res- 
olution of J?~ 110,000. 

We reduced the raw data with the standard IRAF pro- 
cedure and extracted 1-dimensional (ID) spectra. With 
a 20- minutes exposure, we typically achieved a signal-to- 
noise (S/N) of 80 — 90 per pixel of the ID spectra. A high 
S/N template spectrum (S/N- 150) of KOI-94 was ob- 
tained without the Iodine cell for the RV analysis. This 
template spectrum was then subject to deconvolution of 
the instrumental profile (IP), which was reproduced by 
taking the flat lamp spectrum transmitted through the 
Iodine cell. Using this deconvolved template represent- 
ing the intrinsic stellar spectrum, we extracted the RV for 
ea ch of the sp e ctra b y the RV analysis routine developed 
bv lSato et at] (|2002f ). The resulting relative RV's (after 
corrected for the motion of the Earth), along with their 
uncertainties are summarized in Tabled The typical RV 
uncertainty is —10 m s . 

3. ANALYSIS AND RESULTS 
3.1. Stellar Parameters 

Since KOI-94 is not a confirmed planetary system as 
of 2012 September and no spectroscopic measurement 
has been published, we ourselves estimate the basic stel- 
lar parameters from the template spectrum of KOI-94 
used for the RV template. Following lTakeda et al.l ((20021 
120051 ). we estimate the stellar atmospheric parameters 
(T G ff, the surface gravity logg, metalicity [Fe/H], and 
microturbulence dispersion £), based on the equivalent 
widths of Fe I and Fe II lines located at 5000 - 7600 A 
wavelength region. The result of the measurement is 
shown in Table [3] (A). 

We also extract the projected rotational velocity 
UsinJ s by modeling the line profiles of KOI-94's tem- 
plate. By convolving theoretical line profiles for the same 
type star with the rotation plus macroturbulence broad- 
ening kernel and the IP, we fit the line profiles around 
6100 A. In doing so, we assume that the macroturbulence 
dispersion £ is 4.5 ± 1.1 km s , based on the empirical 



TABLE 2 

Radial velocities measured with Subaru/HDS. An 
arbitrary rv offset is subtracted in the data. 



Time [BJD (TDB)] 


Relative RV [m s~ x ] 


Error [m s~ L ] 


24561 09 86607 


-16.9 


14.1 


24561 1 05250 


-14.0 


16.1 


24561 1 84846 


-13.4 


13.1 


24561 1 1 09425 


-24.2 


10.7 


24561 1 1 871 51 


-35.5 


11.1 


24561 1 2 0831 9 


-8.2 


10.9 


24561 1 2 8531 7 


-32.7 


12.7 


24561 1 2 89785 


-31.4 


12.1 


24561 47 97040 


. 8 


9.8 


24561 48 96781 


9 . 7 


13.4 


2456149 75308 


-5.2 


9.6 


2456149 76802 


-13.6 


9.5 


2456149 78275 


4.1 


8.9 


2456149 79747 


14.2 


8.7 


24561 49 81 220 


0.5 


9.6 


24561 49 82692 


15.1 


10.1 


24561 49 841 63 


12.5 


9.6 


24561 49 85635 


16.2 


9.6 


24561 49 871 07 


13.5 


9.7 


24561 49 88579 


9.4 


9.6 


24561 49 90050 


9.9 


8.8 


24561 49 91 522 


-7.1 


9.8 


2456149.92994 


-13.6 


9. 1 


2456149.94466 


-16.4 


10.0 


2456149.96175 


-20.8 


9.3 


2456149.97647 


-34.6 


8.6 


2456149.99119 


-45.8 


9.1 


2456150.00591 


-44.6 


9.2 


2456150.02063 


-43.2 


9.2 


2456150.03536 


-44.0 


9.2 


2456150.05008 


-28.3 


9.8 


2456150.06480 


9.9 


10.4 


2456150.07952 


-3.6 


10.1 



relation bv lValenti fe Fischer! (120051 ) . The resulting best- 
fit value and its uncertainty for V sin I s are also shown in 
Table[3](A). This procedure and validity of our measure- 
ment of V sin I s are described by iHirano et~aT1 ((2012T) 111 
detail. 

Given the atmospheric parameters, we estimate the 
mass and age of the host st ar. We here ad opt the Yonsei- 
Yale (Y 2 ) isochrone model (|Yi et al.ll2001|) . Based on this 
model, we compute the best-fit values and their uncer- 
tainties for the mass M s , radius R s , and age of KOI-94. 
Table [3] (A) also shows the result of those estimates. 

3.2. Modeling and Fitting of the Observed RV's 

The upper panel of Figure Q] plots the RV data phase- 
folded with the orbital period of KOI-94.01 (i.e., ~ 22.34 
days), and the lower panel is the zoomed- in version 
around the transit phase after subtracting the orbit of 
KOI-94.01. The RV variation during the transit of KOI- 
94.01 clearly indicates that the RM effect is securely de- 
tected, and that the planetary orbit is prograde with 
respect to th e host star's e q uator. 

Foll owing IHirano et al.l (|2011al) and iNarita et al.l 
we fit the observed RV's assuming that the RV's 
can be approximated by the RV variation due to Ke- 
plerian motion of KOI-94.01 and the velocity anomaly 
due to the RM effect. Because of the lack of RV data 
points, we here neglect the impacts of KOI-94. 02, 94.03, 
94.04 on RV variations. Based on the transit depths 
of planetary candates reported by Kepler and using a 
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TABLE 3 
System Parameters 



Orbital Phase 




-0.01 



-0.005 

Orbital Phase 



0.005 



0.01 



Fig. 1— (Upper) Phase-folded RV data of KOI-94, ob- 
tained with Subaru/HDS. The best-fit model curve is shown in 
red. (Lower) The RV variation around the planetary transit of 
KOI-94. 01 (blue) and its best-fit model (red solid line). The Kep- 
lcrian orbit is subtracted from the RV's. 



simple relation between the planet mass and radius, 
M p /M e = (i? p /i? ffi ) 2 06 , where M m and R ffi are respec - 
tively the Earth mass and radius (jLissauer et al.ll2011l ). 
we roughly estimate the planet mass and compute the 
RV semi-amplitude exerted by each candidate; the ratios 
of RV semi-amplitudes for KOI-94.02, 03, and 04 to the 
semi-amplitude for KOI-94. 01 are estimated to be 0.17, 
0.25, and 0.04, respectively. Considering that the max- 
imum RV variation in Figure [Tj is approximately 30 m 
s _1 and our error for each RV is no less than 10 m s" 1 , 
one-planet approximation for the RV variation of KOI- 
94 is not a bad approximation. Precise measurements of 



Parameter 

(A) Spectroscopic Parameters 

% H (K) 

log g (dex) 

[Fc/H] (dex) 

£ (km s" 1 ) 

Vsin/ S (km s —1 ) 

M a (M ) 

Rs (Rq) 

age (Gyr) 



Value 



6116 ± 30 
4.123 ±0.055 
-0.01 ±0.04 
1.29 ±0.20 
7.33 ±0.32 

i 9 c-+0.03 
± ' zo -0.04 

1 61+° 11 
1 - D1 -0.12 
o q+0.3 
J -0.2 



(B) Orbital and RM Parameters 
K (m s" 1 ) 
Vsinis (km s _1 ) 

a n 

7offset (m s" 1 ) 

x 2 


14.3 ±4.9 

8.011°;™ 

"611? 
_q q+ 2 - 5 

0.80 


(C) Double Transit Parameters 




T c (1) (BJD) 


2455211.51363 ± 0.00024 


T c (3) (BJD) 


2455211.51790 ±0.00062 


«l 


0.10 ±0.06 


U2 


0.61 ±0.08 


8 (°) 


-1.15 ±0.55 



the impacts KOI-94.02, 94.03, and 94.04 are beyond the 
scope of this letter and will be discussed together with 
further follow-up observations. 

The RV variation including a transit phase is usually 
modeled as 

Vmodd = K[cos(f + w) + ecos(tu)] + Aurm, (1) 

where K, e, and vj are the RV semi-amplitude, orbital 
eccentricity, and argument of periastron, respectively. 
Again, because of the lack of RV data points and the fact 
that KOI-94 is a densely packed multiple planetary sys- 
tem with two giant planets being within ~ 0.3 AU from 
the host star, we decided to fix the orbital eccentricity 
as zero, and put a rough constraint on the RV semi- 
amplitude for KOI-94. 01. Future long-term monitoring 
will the possibility of non-zero eccentricity, improve the 
analysis with non- vanishing eccentricity. 

In order to model the RM vel ocity anomaly A^rm , 
we employ the analytic formula bv lHirano et al.1 (|2011bh . 
The analytic formula is derived for the case that RV's are 
extracted by fitting the stellar template spectrum taken 
outside of the transit to the distorted spectrum during 
a transit, and therefore is more suitable for our calibra- 
tion of the RM velocity anomaly than the classical for- 
mula derived by taking the i ntensity- weighted center of 
the distorted line profile (e.g., lOhta et al.ll2005l ). For the 
stellar parameters required for the RM formula, we adopt 
the macroturbulence dispersion of £ = 4.5 km s _1 , and 
intrinsic Gaussian and Lorentzian line widths of f3 = 2.4 
km s _1 and 7 — 1.0 k m s~\ respec t ively, b ased on the 
empirical estimates by iHirano et al.l (|2011b[ ). As for the 
limb-darkening, we employ the quadratic limb-darkening 
law wi t h ui = 0.40 and 112 = 0.31 based on the table by 
IClaretl (|2004D . Also, we here adopt R p /R s = 0.06856, 
a/R s = 26.1, and I a = 89.33° for the transit parame- 
ters of KOI-94. 01, which are the values provided by the 
Kepler team. 

We fit the observed RV's using Markov Chain Monte 
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Carlo (MCMC) algorithm. The x 2 statistics in our case 
is 

E(^modcl,i — K}bs,i) 2 
^2 ' 



X 



(2) 



where V mo dci,i and V b s ,i ar e the i— th modeled and ob- 
served RV's, and a% is its error. The remaining fitting 
parameters are K, the spin-orbit angle A, Tasini,, and 
the RV offset (zero-point) 7 ff sc t of our dataset listed in 
Table [2j The resulting best-fit values and la errors are 
summarized in Table [3] (B). The red curves in Figured] 
indicate the best-fit models for the whole orbit (upper) 
and the RM effect (lower), respectively. 

The best-fit value of V su\I s purely estimated by the 
RM velocity amplitude (Vsinl s = 8.01±o.73 km s _1 ) 
agrees well with the spectroscopically measured value 
(V s'ml s = 7.33 ± 0.32 km s _1 ), validating our measure- 
ment and modeling of the RM effect. One concern in 
our analysis is that since we neglected the impacts of 
the other three planetary candidates, we might have ob- 
tained a systematically biased estimate for the RV semi- 
amplitude K of KCT94.0f . A spurious RV baseline dur- 
ing the transit phase leads to a systematically biased es- 
timate for the spin-orbit angle A, and sometimes makes 
an aligned system look misaligned, or vice versa. Thus, 
we perform the following test; we remove the RV data 
on the transit night, and fit the remaining RV's with K 
and 7offset being free. Then, using the best-fit value of 
K, we fit the RV's during the transit night in order to 
derive A and Vsin/ S for the fixed K. This treatment 
can remove possible systematics caused by instrumental 
effects, the impact of other planets, and/or non-zero ec- 
centricity. As a result, we find Vsin/ S = 8.O2+Q 74 km 

s _1 and A = — 10+ n deg, implying a good spin-orbit 
alignment again. Based on this test, we conclude that 
the measurement of the spin-orbit angle A is robust and 
the orbital axis of KOI-94.01 is almost parallel to the 
projected stellar spin axis. 

Substituting the best-fit value of K and our estimate 
for the stellar mass M s , we obtain ~ 0.23M,/ for the 
planet mass of KCT94.0f. This mass is slightly small 
compared to the theoretical expectation for a plane t hav- 
ing a radius of ~ 9.25i? e (|Mordasini et al.ll2012D . but 
is rather similar to low-density planets like Kcpler-9c 
(|Holman et al.ll2"010D . In our analysis, however, we com- 
pletely neglected the impacts of the other three planet 
candidates so that this estimate of the planet mass more 
or less has systematics. Again, a long-term RV moni- 
toring is required in order to put tight constraints on 
masses and eccentricities of all the planet candidates in 
this system. 

4. DISCUSSION AND SUMMARY 

While we have seen that KOI-94.0I's orbital axis is 
aligned with the stellar rotation axis, it is not evident 
that the orbits of the other planet candidates are also 
neatly aligned with the stellar spin. Fortunately, how- 
ever, we found a lucky lightcurve in which a "double tran- 
sit" event is observed in Kepler's public lightcurve. The 
upper panel of Figure [5] shows the normalized lightcurve 
of KOI-94 around BJD = 2455211.5. In this figure, the 
first transit event starts around BJD = 2455211.35, and 
after a while, it shows the second, deeper dimming. Ac- 
cording to the public Kepler planet candidate ephemeris, 




Fig. 2. — (Upper) The public lightcurve taken by Kepler around 
BJD = 2455211.5. In addition to the simultaneous double transit 
event, there is a bump around the transit center of KOI-94.01, 
which most likely represents the planet-planet eclipse. (Lower) 
The definition of the mutual inclination S between KOI-94.01 and 
94.03. The orbits of the two planets are projected onto the sky 
plane. This is just a schematic description and does not reflect the 
the real planet sizes and angle between the planetary orbits. 

the first shallow dimming represents the transit of KOI- 
94. 03 and the second deeper one corresponds to a tran- 
sit of KOI-94.01. Interestingly enough, there is a bump 
around the transit center with height of ~ 0.0003 in 
terms of the relative flux. Similar events have been fre- 
quently seen when the trans iting planet crosses a starspot 
(|Sanchis-Oieda et~ai1 1201 lh . We checked the out-of- 
transit lightcurve (long-cadence data) and computed the 
periodogram in order to look for a signal induced by 
starspots. Consequently, we found a weak peak around 
the period of ~ 11 days. Howev er, the flux variabil ity 
in the lightcurve (as is defined in iHirano et al.l 1201 2h is 
about 0.02% so that the size of the bump seen in Figure 
[2] cannot be explained by a spot-crossing event. Consid- 
ering the timing and size of the bump, the most likely 
scenario would be the eclipse of the planet by another 
transiting pla net (planet-planet eclipse, al so called "mu- 
tual event", iRagozzine k Holmanl 120101) . Since KOI- 
94. 03 has the longer orbital distance from the host star, 
a part of the inner planet, KOI-94.01, is occulted by 
the outer transiting planet, KOI-94. 03, during the eclipse 
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event represented by the bump. 

Since the impact parameters of the two transiting plan- 
ets are fixed from the transit photometries, the only pa- 
rameter to determine the relative planetary orbits on the 
stellar disk and describe the bump in Figure [2] is the mu- 
tual inclination between the two transiting planets. If 
the mutual inclination of the two planets is considerably 
large, we expect no planet-planet eclipse. The timing and 
shape of the anomalous bump in Figure [2] put a strong 
constraint on the mutual inclination. 

We define the mutual inclination projected onto the 
sky as S (sec the lower panel of Figure [2|) and model the 
anomalous bump by a simple geometric calculation as a 
function of time and 5. The mutual inclination S sensi- 
tively changes the duration of the planet-planet eclipse, 
and can be constrained by the shape of the bump around 
the bottom of the transit lightcurve; the resulting bump 
is longer for nearly prograde [6 ~ 0), and shorter for ret- 
rograde planets (5 « ir). By modeling the double transit 
event as a sum of two single transit lightcurves and the 
bump function (which is expressed analytically as an area 
fraction of the overlapping region by the two transiting 
planets), we fit the observed lightcurve by MCMC algo- 
rithm. The free parameters in our fit are the mid-transit 

times T c (1) and T c (3) for KOI-94.01 and KOI-94.03, the 
limb-darkening parameters m and ui for the quadratic 
limb-darkening law 6 , and the projected mutual inclina- 
tion <5. The other transit parameters (i.e., R P /R S7 a/R s , 
and i ) for both of 94.01 and 94.03 are fixed at the values 
delivered by the Kepler team assuming circular orbits for 
both of the candidates. Table [3] (C) shows the best-fit 
values and their uncertainties of our fitting. The red solid 
line in the upper panel of Figure [2] indicates the best-fit 
model and the residual of observed data from the best-fit 
model is shown at the bottom. The orbital planes of the 
two transiting planets are remarkably well aligned, which 
means that host star's spin axis and planetary axes of 
both KOI-94.01 and 94.03 are all well aligned. 

Like the Keplcr-30 system, the spin-orbit alignment 
in the KOI-94 system implies that those planets in this 
multiple transiting system have experienced a quiescent 
migration process rather than having been carried to 
the present locations by dynamical scattering processes 
or long term perturbations by outer objects. This has 
also ruled out the possibility that an earlier magnetic 
interaction between the star and disk causes a grad- 
ual misalignmen t of the stellar spin from the disk plane 
(|Lai et al.l 120 1 If) at least for this system. On the other 



hand, our result should reinforce the hypothesis that 



spin-orbit misalignments are only seen for "isolated" hot 
Jupiters, which may have m igrated by chaotic processes 
(jSanchis-Oieda et all 120121) . 

Finally, we note that KOI-94 has the effective temper- 
ature of T ff = 6116 ± 30 K, which suggests that the star 
most likely has a convective envelope rather than being 
fully covered by a thi ck radiative envelope as seen for 
massive (hotter) stars (|Pinsonneault et all 1200 1|) . Thus, 
this system does not provi de us the chance to test the 
recent theoretical work by iRogers et al.l (|2012h . stating 
that the stellar spins of massive stars with radiative en- 
velopes and convective cores may spontaneously change 
directions due to the internal gravity waves (IGW) gen- 
erated around the border between the radiative envelope 
and convective core, regardless the presence of planets. A 
future observation of the RM effect for a multiple tran- 
siting system with a massive host star will confirm or 
refute this intriguing scenario. 
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